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串联 回转 体高 速 和 人 水 空 化 及 运动 特性 的 数值 模拟 
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摘 要 :为 了 研究 串联 回转 体高 速 入 水 问题 ,基于 求解 雷诺 时 均 的 Navier-Stokes 方程 ,采用 VOF 多 
相 流 模型 和 6 自由 度 (6DOF) 刚 体 运 动 模型 ,引入 Schnerr-Sauer 空 化 模型 ,结合 重 登 网 格 技术 建立 
了 串联 回转 体高 速 入 水 的 三 维 数值 计算 模型 ,发现 了 串联 回转 体 入 水 空 泡 演 化 的 一 些 新 现象 ,获得 
了 不 同 入 水 时 序 下 ee ds an 结合 入 水 空 泡 的 演化 规律 ,分 析 了 入 水 空 泡 演 

化 过 程 对 串联 回转 体 运动 特性 的 影响 ,给 出 了 串联 回转 体高 速 入 水 的 运动 模式 。 结 果 表 明 : 次 发 回 
转 体 入 水 过 程 受 gg 在 入 水 撞击 阶段 ,次 发 回转 体 受 到 的 冲击 力 小 
于 首发 回转 体 的 ,而 且 在 入 水 过 程 中 次 发 回转 体 的 速度 衰减 小 于 首发 回转 体 的 ; 随 着 入 水 时 序 增 
大 ,串联 回转 体 入 水 运动 模式 依次 出 现 追尾 模式 和 弱 扰 动 模式 ,在 追尾 模式 中 ,次 发 回转 体能 够 进 
入 首发 回转 体形 成 的 入 水 空 泡 内 部 运动 并 最 终 追 上 首发 回转 体 ,在 弱 扰 动 模式 中 ,首发 回转 体 的 入 
水 空 泡 演化 过 程 不 受 次 发 回转 体 入 水 的 影响 ,而 次 发 回转 体 的 入 水 空 泡 演 化 和 运动 受到 了 首发 回 
转 体 空 泡 尾 迹 的 干扰 。 
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Numerical simulation of cavitation and motion characteristics 
of high-speed water entry of tandem revolved bodies 
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Abstract :In order to study the high-speed water entry of tandem revolved bodies ,a three-dimensional nu- 
merical simulation model of high-speed water entry of tandem revolved bodies is established based on sol- 
ving Reynolds time averaged Navier Stokes equation ,using VOF multiphase flow model and 6 degrees of 
freedom (6DOF ) rigid body motion model, Schnerr-Sauer cavitation model and overlapping grid technolo- 
gy. Some new phenomena in the cavities evolution of tandem revolved bodies were found ,and the evolution 
law of cavities under different water entering time sequences were obtained. Combined with the evolution 
law of cavities ,the influence mechanism of the evolution process of cavities on the motion of tandem re- 


volved bodies was analyzed ,and the motion mode of tandem revolved bodies entering water at high speed 
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was given. The results show that the water entry process of the secondary revolved body is affected by the 


cavitation flow field of the first revolved body. In the water entry impact stage ,the impact force of the sec- 


ondary revolved body is less than that of the first revolved body. Moreover ,the velocity attenuation of the 


secondary revolved body is less than that of the first revolved body. With the increase of water entering 


time sequence ,the water entry motion mode of tandem revolved bodies appears rear-end collision mode and 


weak disturbance mode in turn. In the rear-end collision mode ,the secondary revolved body can enter the 


water entry cavity formed by the first revolved body and finally catch up with the first revolved body. In the 


weak disturbance model ,the cavities evolution process of the first revolved body is not affected by the wa- 


ter entry of the second revolved body ,while the cavities evolution and motion of the second revolved body 


are disturbed by the cavitation wake of the first revolved body. 


Key words :tandem revolved bodies; high speed water entry; cavity; numerical simulation 


在 实际 军事 应 用 中 ,为 了 提高 武器 的 打击 火力 、 
命中 效率 和 毁伤 效果 ,往往 需要 在 短 时 间 内 连续 发 
射 超 空 泡 射 弹 / 反 潜 火 箭 进行 攻击 。 例 如 ,水 下 蛙 人 
或 天 下 机 器 人 的 移动 速度 通常 为 0.5 ~5 m/s,10 ms 
内 这 些 目标 的 位 移 也 只 有 几 厘 米 " ,那么 朝 同 一 个 
六 而 连续 发 射 射 弹 可 以 提高 目标 被 消灭 的 概率 。 上 
述 凶 程 应 用 涉及 了 串联 运动 体 的 高 速 人 水 问题 ,该 
人 炒 过 程 不 仅 包含 穿越 自由 面 滑动 `. 相 变 等 复杂 流 
动 状态 ,还 存在 各 运动 体 流 场 之 间 的 相互 干扰 ,水 动 
和 | 性 十 分 复杂 ,这 与 单独 的 运动 体高 速 入 水 有 很 
大 多 区 别 。 因此 ,研究 串联 运动 体高 速 和 人 水 问题 具 
存 刘 要 的 科学 意义 和 工程 应 用 价值 。 

会 目前 国内 外 学 者 对 入 水 问题 的 研究 主要 集中 在 
单产 物体 和 水” 。1897 年 , Worthington 等 中 采用 内 
光 鳗 影 技 术 首次 获得 了 球体 低速 人 水 产生 的 物理 现 
ee A er 


验 ， 二 s 泡 产生 的 影响 以 
及 射 弹 运 动 参数 对 入 水 空 泡 发 展 的 影响 。Abel- 

ont 中 对 射 弹 以 不 同 角度 入 水 (入 水 角度 分 别 为 45°、 
60° .90°) 进行 了 实验 研究 ,测量 了 入 水 空 泡 的 压力 。 
Lee 等 根据 能 量 守 恒 原 理 建立 了 模拟 射 弹 高 速 入 
水 过 程 中 空 泡 发 展 和 演 灭 的 理论 模型 ,可 以 预测 入 
水 空 泡 深 闭合 的 时 间 和 位 置 。Shi 等 中 采用 高 速 报 
像 技 术 记 录 了 射 弹 高 速 垂直 和 人 水 的 图 像 , 获 得 了 射 
弹 入 水 后 速度 的 变化 ,并 提出 了 一 个 理论 模型 来 预 
测 射 弹 入 水 后 的 运动 。Truscott 等 1 进行 了 射 弹 小 
角度 倾斜 人 水 实验 ,为 了 使 射 弹 能 够 成 功 人 水 并 在 
水 下 稳定 运动 ,他 们 重新 设计 了 射 弹头 部 的 几何 结 
构 ,讨论 了 弹头 的 几何 形状 对 射 弹 的 动力 学 .临界 人 
射 角 、 入 水 空 泡 发 展 的 影响 。Iranmanesh 等 利用 


VOF 方法 对 水 平 圆 柱 入 水 的 流动 特性 进行 了 三 维 数 
值 模拟 ,研究 了 圆柱 体 的 直径 、 长 度 、 入 水 速度 等 参 
数 对 入 水 深度 的 影响 。Nguyen 等 (中 看 合 Navier- 
Stokes( N-S) 方 程 和 6 自由 度 (6DOF) 刚 体 运动 模 
型 ,基于 重 辣 网 格 技 术 建 立 了 求解 物体 入 水 的 三 维 
数值 模拟 方法 ,并 验证 了 该 方法 的 可 行 性 。Treichler 
等 "采用 高 速 摄像 技术 实验 研究 了 超 空 泡 射 弹 垂 
直人 水 引起 的 空 化 流动 特性 ,并 采用 压 电 传感器 记 
录 了 撞击 点 位 置 和 射 弹 运动 方向 上 的 压力 。Yao 
等 '5] 通 过 实验 和 理论 相 结合 的 方法 研究 了 弹丸 垂 
直入 水 问题 ,获得 了 入 水 空 泡 的 发 展 过 程 和 弹丸 的 
运动 特性 ,并 指出 射 弹 在 撞 水 之 后 的 阶段 需要 考虑 
重力 效应 。Jiang 等 '" 通过 实验 和 数值 模拟 的 方法 
研究 了 射 弹 垂直 进入 减 阻 液 中 诱导 的 超 空 泡 流 动 特 
性 。 马 庆 鹏 等 "采用 VOF 多 相 流 模型 建立 了 锥 头 
圆柱 体高 速 和 人 水 的 二 维 数值 模型 , 获得 了 运动 体 的 
入 水 速度 变化 及 入 水 空 泡 的 生成 与 发 展 过 程 。 张 伟 
等 "利用 轻 气 炮 进行 了 平头 和 球形 两 种 弹 体 水 平 
人 水 实验 ,获得 了 弹 体 和 人 水 产生 的 冲击 波 和 传播 特 
性 。 时 素 果 等 ”在 水 池 中 进行 了 细 长 体高 速 入 水 
实验 ,获得 了 入 水 空 泡 的 形态 和 细 长 体 的 表面 压力 
变化 规律 。 郭 子 涛 等 …” 利用 一 级 轻 气 炮 设备 开展 
了 多 种 头 型 的 弹 体 水 平和 水 实验 ,入 水 速度 范围 为 
100 ~ 300 m/s, 获 得 了 入 水 空 泡 的 演化 过 程 和 空 泡 
壁 的 扩展 速度 特性 。Chen 等 “| 实验 研究 了 射 弹 高 
速 入 水 过 程 的 弹道 稳定 性 ,分 析 了 弹头 形状 、 撞 击 速 
度 、 入 水 姿态 角 对 射 弹 弹 道 稳 定性 的 影响 ,结果 显示 
平头 射 弹 的 弹道 稳定 性 最 好 。 王 晓 辉 等 ”| 基于 灶 
合 欧 拉 - 拉 格 朗 日 方法 开展 了 射 弹 高 速 人 水 过 程 的 
数值 模拟 ,分 析 了 射 弹 入 水 尾 拍 这 一 动态 运动 的 形 
成 机 理 。 
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关于 串联 运动 体 人 水 的 研究 目前 尚 不 多 见 。 何 
春 涛 等 ”对 在 重力 加 速 下 两 个 弹 体 串联 倾斜 低速 
人 水 问题 进行 了 探索 试验 ,结果 显示 两 弹 体 间隔 距 
离 较 小 的 情况 下 发 生 了 追尾 碰撞 ,碰撞 导致 两 弹 体 
偏离 了 预期 运行 轨迹 。 余 德 硕 等 .” 利用 高 速 摄像 
技术 进行 了 串联 回转 体 低 速 垂 直入 水 实验 ,得 到 了 
深 闭合 空 泡 扰动 下 串联 入水 的 运动 模式 。Yun 等 '" 
采用 高 速 摄像 技术 开展 了 串联 双 球 体 低 速 倾斜 人 水 
实验 ,获得 了 入 水 空 泡 的 演化 特性 、 球 体 的 运动 特性 
及 受 力 变化 。 施 红 辉 等 开展 了 两 连 发 超 空 泡 射 
弹 垂 直人 水 的 二 维 数值 模拟 , 发 现 了 第 二 发 射 弹 的 
流 场 会 受到 第 一 发 射 弹 入 水 流 场 的 干扰 。 

综 上 所 述 ,国内 外 对 单独 物体 入 水 研究 较为 成 
熟 , 但 关于 串联 物体 高 速 入 水 问题 的 研究 仍 处 于 起 
步 阶段 。 串 联 物体 高 速 入 水 的 流 场 非常 复杂 ， 
了 T 吨 介 质 滑动 ` 空 化 ` 流 场 的 相互 作用 ,其 人 水 空 
的 演化 机 理 流体 动力 学 特性 等 尚 不 明确 , 还 需 
-上 的 深入 研究 。 本 研究 基于 求解 N-S 方程 ,采用 
V 怠 多 相 流 模型 和 6DOF 刚体 运动 模型 , 引入 
Sehuerr- Sauer 空 化 模型 ,结合 重 麦 网 格 技术 实现 流 
场 和 弹道 的 看 合 求解 ， 建立 串联 回转 体 垂 直 高 速 人 
大 六 数值 模型 ,获得 不 同人 水 时 序 下 串联 回转 体 入 
水 宝 泡 的 演化 规律 和 运动 特性 ,分 析 新 的 流动 现象 
ne 


控制 方程 与 数值 计算 方法 


1 无 控制 方程 


二 本 研究 基于 均 质 混合 流 假设 处 理 和 水 过 程 的 多 
相 流 场 , 即 空气 \ 水 \ 水 蒸气 三 相 的 混合 流体 被 当 作 
密度 和 黏度 随时 间 和 空间 变化 的 一 种 均 质 流体 。 采 
用 流体 体积 函数 (VOF ) 多 相 流 模 型 来 捕获 不 同 流体 
相 之 间 的 界面 。 连 续 性 方程 和 动量 方程 分 别 为 

dp» 9(Cpati) 
et 
其 中 :pw 为 混合 相 密度 ,p。= ap + Qsp, + oipi;a 
as ai 分 别 为 水 蒸气 .空气 ,水 的 体积 分 数 ,a, + a + 
oa =1,p, spl 分 别 为 水 芋 气 .空气 水 的 密度 ; 1 为 
时 间 ;w; 为 i 轴 方 向 的 速度 。 


=0 (1) 


0( pu;) oPnuit) = 
ot Ox; 
oP 0 Ou. ou; 
一 二 一 二 十 2 
i jo sl | pngi (2) 
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其 中 :P 为 压力 ;ww 为 混合 相 动 力 黏度 ,ps = ay + 
Qs + ay Was Wi 分 别 为 水 蒸气 、 空 气 , 水 的 动力 
黏度 ;iu 为 j 轴 方向 的 速度 ; g; 为 重力 加 速度 的 
分 量 。 


采用 SST k-w 消 流 模型 “对 RANS 方程 提供 漠 
流 封闭 ,该 模型 考虑 了 满 流 剪 切 应 力 传输 。 运 动 体 
高 速 人 水 过 程 中 会 产生 空 化 现象 ,对 于 空 化 的 求解 
采用 Schnerr-Sauer 空 化 模型 ,该 空 化 模型 具有 较 
高 收敛 速度 和 计算 稳定 性 。 其 控制 方程 如 下 。 


(op) +Y， opsv,) =R -RR (3) 


PP EN 

Re | = | Pp, >P (4) 

= 3 /2/2-P, 

Ree tl | I <P (5) 
其 中 :为 水 蒸汽 相 的 速度 矢量 ;R。R. 分 别 为 蒸发 
速率 和 冷凝 速率 ; R, 为 气 核 的 半径 , R = 

Q, 3 1 1 要 11 。 过 人 藻 
放下 村 省 =1x10" ;为 水 的 饱和 大 
汽 压力 。 


物体 的 运动 通过 6DOF 刚体 运动 模型 来 求解 。 
物体 运动 可 分 为 质心 的 平 动 和 绕 质 心 的 转动 ,在 每 
个 时 间 步 ,利用 作用 在 物体 的 力 和 力矩 计算 出 平移 
加 速度 和 角 加 速度 ,然后 通过 积分 得 到 物体 的 平移 
速度 和 角速度 ,进一步 计算 得 到 物体 新 的 质心 位 置 
和 欧 拉 角 '“ 。 设 六 和 妨 .分 别 表示 第 二 时 间 步 质心 
的 位 置 和 物体 的 运动 方向 ,那么 第 n+1 时 间 步 质心 
的 位 置 和 物体 的 运动 方向 为 
XY = + At 
I =0", +Gw.At 
其 中 :At 为 时 间 步 长 ;w 和», 分别 为 物体 的 角速度 
和 平移 速度 ; G 为 变换 矩阵 。 


1.2 数值 模型 及 边界 条 件 
回转 体 模型 采用 参考 文献 [19 ] 中 的 平头 弹 体 ， 


如 图 1 所 示 。 弹 体 的 直径 D =6 mm, 长 度 工 = 
24 mm ,材质 为 45# 钢 , 质 量 为 4. 88 g。 


L=24 mm 


(6) 


6 mm 


D= 


妈 1 回转 体 结构 示意 图 


Fig.1 Schematic diagram of the revolved body model 
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计算 域 及 边界 条 件 设置 如 图 2 所 示 。 计 算 域 尺 
度 为 50D x 50D x250D, 流 场 空 气 域 高 度 为 25D, 水 
域 深度 为 225D, 首 发 回转 体 与 次 发 回转 体 首 尾 相 中 
12D, 计 算 域 上 方 为 压力 入 口 ,压强 与 标准 大 气压 一 
致 ,计算 域 侧 面 和 底部 为 压力 出 口 ,它们 的 压力 值 通 
过 自 定义 函数 进行 指定 ,首发 回转 体 和 次 发 回转 体 
的 表面 均 设 定 为 无 滑 移 壁面 。 


25D 


225D 


压力 出 口 


图 2 计算 域 及 边界 条 件 设置 


Fi2 Computational domain and boundary condition setting 
¢ 


1.3= 网 格 划分 及 数值 方法 

GO 流 场 计算 采用 重 琶 (overset) 网 格 技术 ,图 3 给 
出 了 计算 域 的 网 格 划分 示意 图 。 该 计算 域 网 格 划 分 
分 为 背景 网 格 (background mesh) 和 前 景 网 格 (com- 
ponent mesh) ,由 于 首发 回转 体 和 次 发 回转 体 之 间 存 
在 相对 运动 ,其 前 景 网 格 计算 域 应 当 相互 独立 ,因此 
需要 划分 两 套 前 景 网 格 。 图 3 中 3 套 网 格 均 采 用 六 
面体 结构 化 网 格 , 网 格 数量 约 为 214 万 。 背 景 网 格 
区 域 对 自由 液 面 附近 和 回转 体 运动 区 域 进行 局 部 加 
密 , 提 高 人 水 空 泡 边 界 的 计算 精度 。 两 套 前 景 网 格 
区 域 的 大 小 和 网 格 划分 方式 相同 , 均 为 包 庄 回 转 体 
的 圆柱 ,其 长 度 为 10D, 直 径 为 5D, 前 景 网 格 区 域内 
部 的 回转 体 壁 面 附件 的 网 格 进行 局 部 加 密 。 计 算 中 ， 
设置 前 景 网 格 1 计算 域 与 首发 回转 体 同步 运动 以 及 
前 景 网 格 2 计算 域 与 次 发 回转 体 同步 运动 ,然后 结合 
6DOF 算法 可 以 实现 两 个 回转 体 各 自 独立 运动 的 六 自 
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由 度 运 动 求解 。 


前 景 网 格 1 
图 3 网 格 划 分 示意 图 
Fig.3 Schematic diagram of computational domain grids 

数值 模拟 基于 Flunet 19.0 平台 开展 ,其 中 瞬时 
压力 场 与 速度 场 的 耦合 求解 采用 Coupled 算法 , 压 
力 插 值 采用 PRESTO! 格式 ,各 项 体积 率 离散 采用 
Modified HRIC 格式 , 湛 流 和 耗 散 相 离 散 均 采用 二 阶 
迎风 格式 。 

为 排除 网 格 密度 对 数值 模拟 结果 的 影响 ,需要 
开展 网 格 无 关 性 验证 。 在 自由 液 面 附近 、 回 转 体 壁 
面 和 运动 区 域 进行 了 局 部 加 密 ,建立 了 3 种 不 同 密 
度 的 网 格 ,网 格 数量 分 别 约 为 102 万 (casel) 、214 万 
(case2) 307 万 (case3 ) 。 对 和 水 时 序 为 1 ms, 和 人 水 
速度 为 106. 8 m/s 的 串联 回转 体 入 水 进行 了 数值 模 
拟 , 得 到 了 首发 回转 体 运 动 速度 随时 间 的 变化 曲线 ， 
如 图 4 所 示 。 从 图 中 可 以 看 出 ,casel 与 case2 网 格 
计算 结果 差异 较 大 ,而 case2 和 case3 网 格 的 计算 结 
果 基 本 一 致 。 综 合 考虑 计算 效率 和 计算 精度 ,计算 
时 选择 了 case2 的 网 格 开展 数值 计算 。 
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图 4 不 同 网 格 密度 的 计算 结果 


Fig.4 Calculated results of different mesh densities 
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1.4 数值 方法 验证 


关于 串联 物体 高 速 入 水 目前 还 未 见 到 相关 的 实 
验 数据 ,因此 采用 上 文 的 数值 模拟 方法 对 文献 [19] 
中 实验 开展 的 单 发 平头 弹 体 垂直 和 人 水 问题 进行 数值 
模拟 。 本 研究 选择 的 回转 体 模型 ( 见 图 1) 就 是 该 实 
验 中 所 用 的 平头 弹 体 模型 ,实验 中 平头 弹 体 的 人 水 
初速 度 为 106. 8 m/s。 数 值 模 拟 计算 过 程 中 ,设置 次 
发 回转 体 静止 不 动 ,首发 回转 体 以 106. 8 m/s 的 速 
度 入 水 ,这 样 就 可 以 实现 单独 的 物体 入 水 问题 求解 。 
站 从 出 了 回转 体 人 水 空 泡 演化 的 实验 和 数值 模拟 

结果 。 通 过 对 比 可 以 看 出 两 者 对 应 良好 ,单独 回转 
体 入 水 经 历 了 入 水 撞击 ` 开 空 泡 、 空 泡 表 面 闭合 、 空 
泡 随 动 及 尾部 溃 灭 等 阶段 。 在 实验 和 数值 模拟 中 均 
表明 该 回转 体 人 水 的 弹道 稳定 性 很 好 ,几乎 没有 发 
红叶 道 的 偏 移 和 弹 体 姿 态 角 的 变化 。 图 6 给 出 了 加 
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Fig.$5 Comparisons of evolution of water-entry cavities 


between numerical simulations and experiments 


用 力学 学 


转 体 入 水 深度 随时 间 的 变化 ,将 数值 模拟 结果 和 实 
验 结果 进行 对 比 , 发 现 两 者 具有 较 好 的 一 致 性 。 通 
过 上 述 对 比 ,说 明 该 数值 模拟 方法 是 有 效 的 。 


-@ 数值 模拟 结果 


让 SO 3.00 3:3 “0 
tms 


到 6 单独 回转 体 入 水 深度 随时 间 的 变化 


Fig.6 Variation of water entry depth 


of single revolution body with time 
2 结果 与 讨论 


将 首发 回转 体 头 部 刚 接触 静止 自由 液 面 的 时 刻 
设 定 为 时 间 /的 零点 。 首 发 回转 体 的 人 水 初速 度 为 
Jo ,次 发 回转 体 的 人 水 初速 度 为 Po ,首发 和 次 发 回 
We 自由 液 面 的 时 间 差 定义 为 人 水 时 序 

。 数 值 模拟 中 可 以 通过 UDF 控制 次 发 回转 体 在 
| 
的 入 水 时 序 入 水 ,入 水 过 程 中 ,Vio = Vo =106.8 m/s。 


2.1 典型 工 况 中 串联 回转 体 入 水 空 泡 形态 演化 
特性 


图 7 给 出 了 Ar=lms 时 ,串联 回转 体 和 人 水 空 泡 
形成 与 演化 过 程 。 在 首发 回转 体 开 空 泡 阶 段 , 次 发 
回转 体 紧 随 而 来 ,之 后 次 发 回转 体 进入 到 先 发 射 弹 
的 入 水 空 泡 内 部 ,在 表面 张力 和 水 域 压力 的 作用 下 ， 
入 水 空 泡 在 2 ms 时 刻 发 生 表面 闭合 。 由 于 空 泡 内 
部 主要 为 空 :化 生成 的 区 气相 和 外 界 进入 的 空 气相 ， 
其 密度 远 小 于 水 的 密度 ,那么 次 发 回转 体 受到 的 阻 
力也 远 小 于 首发 回转 体 的 阻力 。 随 着 时 间 的 增加 ， 
两 发 回转 体 之 间 的 距离 越 来 越 小 ,在 2.85 ms 时 ,次 
发 回转 体 追 上 首发 回转 体 将 会 发 生 碰撞 , 这 里 需要 
说 明 的 是 本 研究 没有 研究 两 发 回转 体 碰撞 后 的 运动 
情形 , 磁 撞 涉及 较 复 杂 的 力学 机 制 ,还 需要 进一步 的 
深入 研究 。 
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图 At=1 ms 时 ,串联 回转 体 入 水 空 泡 形成 与 演化 过 程 


Fig.7 When At=1 ms,the formation and evolution 


of water-entry cavities of tandem revolved bodies 
图 8 给 出 了 At=3 ms 时 ,串联 回转 体 人 水 空 

形成 与 演化 过 程 。 首 发 回转 体 的 入 水 空 
面团 合 和 被 拉 离 自由 面 时 ,次 发 回转 体 才刚 接触 自 
插 前 (图 中 := 3 ms 时 ) ,之 后 次 发 回转 体 经 历 了 入 
术 挤 击 \. 开 空 泡 。 随 着 次 发 回转 体 进一步 向 下 运动 ， 

次 发 回 转 体 逐渐 进入 到 首发 回转 体 的 入 水 空 泡 内 
部 533f 引 起 两 个 人 水 空 泡 的 融合 ,融合 后 的 空 泡 在 水 
域 奈 力 和 表面 张力 的 作用 下 在 自 由 面 附近 发 生 了 表 
面 闭合 。 随 着 时 间 的 推移 ,次 发 回转 体 所 在 位 置 周 
辐 的 空 ; 泡 发 生前 注 ， 导致 次 发 回转 体 跟 随 空 泡 的 尾 
部 运动 ,如 图 中 1;=5.4 ms 所 示 。 由 于 次 发 回转 体 的 
速度 高 于 首发 回转 体 ,次 发 回转 体 在 :=7 ms 时 再 次 
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图 8 At=3 ms 时 , 申 联 回转 体 入 水 空 泡 形 成 与 演化 过 程 


Fig.8 When At=3 ms,the formation and evolution 


of water-entry cavities of tandem revolved bodies 
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完全 进入 到 首发 回转 体 的 空 泡 内 部 。 在 1=7.2 ms 
时 ,次 发 回转 体 与 首发 回转 体 发 生 追 尾 。 

图 9 给 出 了 At =5 ms 时 ,串联 回转 体 和 人 水 空 泡 
形成 与 演化 过 程 。 首 发 回转 体 入 水 过 程 先后 经 历 了 
入 水 撞击 、 开 空 泡 .表面 闭合 、 空 泡 随 动 及 演 灭 ,其 空 
泡 演化 不 再 受到 次 发 回转 体 入 水 的 影响 ,与 单 发 回 
转 体 和 人 水 相同 。 而 次 发 回转 体 人 水 过 程 中 受到 了 首 
发 回转 体 空 泡 尾 迹 的 干扰 ,其 入 水 空 泡 形态 与 单 发 
回转 体 和 水 空 泡 存在 明显 区 别 ,而 且 较 早 发 生 颈 缩 
和 省 灭 现象 。 
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图 9 Al =5 ms i 串联 回转 体 入 水 空 泡 形成 与 演化 过 程 


Fig.9 When At=5 ms,the formation and evolution 


of water-entry cavities of tandem revolved bodies 

图 10 表示 的 是 次 发 回转 体 人 水 深度 HH=7D 
时 ,不 同人 水 时 序 下 回转 体 的 压力 云图 。3 种 时 序 
下 ,首发 回转 体 均 在 其 头 部 前 端 产生 高 压 区 ,其 压力 
场 呈 球面 形式 向 周 向 传播 ,低压 区 出 现在 人 水 空 泡 
内 部 ,由 于 空 化 效应 其 压力 值 等 于 水 的 饱和 蒸汽 压 
( 约 为 3540 Pa) 。 入 水 时 序 为 1 ms 时 ,次 发 回转 体 
被 首发 回转 体 的 低压 空 泡 包围 ,其 头 部 前 端 无 高 压 
区 产生 ,因此 没有 形成 人 水 空 泡 。 入 水 时 序 为 3 ms 
时 ,次 发 回转 体 头 部 前 端 产生 了 高 压 区 ,最 大 压力 值 
约 为 15.6 个 标准 大 气压 ,并 形成 了 低压 人 水 空 泡 
但 受到 首发 回转 体 空 泡 尾部 的 影响 ,次 发 回转 体 头 
部 前 方 高 压 区 的 范围 受到 了 压缩 ,不 再 呈 球 面 形式 
向 周 向 传播 。 入 水 时 序 为 5 ms 时 ,首发 回转 体 的 尾 
部 空 泡 发 生 演 灭 并 产生 局 部 高 压 , 在 该 入 水 深度 次 
发 回转 体 受到 首发 回转 体 流 场 的 干扰 较 弱 ,次 发 回 
转 体 正常 形成 了 入 水 空 泡 并 在 其 头 部 前 端 产生 了 呈 
球面 形式 传播 的 高 压 区 ,其 最 大 压力 约 为 33.8 个 标 
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准 大 气压 。 发 回转 体 速度 衰减 程度 大 于 次 发 回转 体 的 速度 衰 


(a) At=1 ms 


(b) A[=3 ms 


二 (c) As ms 

图 10 不 同人 水 时 序 下 ,次 发 回转 体 

OO 入 水 深度 委 =7D 时 的 压力 云图 

QQ 10 Contours of pressure when the water entry depth 
of the secondary revolved body H =7D under 

CN different water entering time sequences 

AN 


2 人 串联 回转 体 入 水 运动 特性 


二 图 11 表示 的 是 入 水 时 序 为 1 ms 时 ,串联 回转 
量 纲 入 水 深度 .速度 .加 速度 随时 间 的 变化 曲 
线 s 自 由 液 面 为 HH=0 的 位 置 。 图 中 瓦 ao 参数 的 
方面 为 竖 直 方向 。 从 图 11(a) 中 可 以 看 出 , 随 着 次 
发 三 转 体 的 入 水 ,两 发 回转 体 之 间 的 距离 越 来 越 小 ， 
在 @ 85 ms 时 ,次 发 回转 体 追 上 首发 回转 体 。 从 图 
11(b) 和 (c) 中 可 以 看 出 ,首发 回转 体 和 人 水 起 始 阶段 
加 速度 较 大 ,速度 衰减 很 快 ,这 是 由 于 物体 撞击 自由 
液 面 时 ,由 于 介质 突变 产生 巨大 的 反 向 冲击 作用 ,使 
得 回转 体 产生 巨大 的 反 向 加 速度 ,速度 衰减 很 大 。 
之 后 随 着 撞击 点 距离 的 增加 和 入 水 空 泡 的 形成 , 首 
发 回转 体 加 速度 值 迅速 降低 。 对 于 次 发 回转 体 , 由 
于 次 发 回转 体 人 水 过 程 中 直接 进入 了 首发 回转 体 的 
入 水 空 泡 内 ,导致 次 发 回转 体 受到 量 级 很 小 的 空气 
阻力 值 ,引起 次 发 回转 体 的 加 速度 和 速度 变化 很 小 。 

图 12 表示 的 是 入 水 时 序 为 3 ms 时 ,串联 回转 
体 无 量 纲 人 水 深度 .速度 .加 速度 随时 间 的 变化 曲 
线 。 从 图 中 可 知 ,次 发 回转 体 在 3 ms 时 刻 人 水 ,入 
水 后 次 发 回转 体 与 首发 回转 体 之 间 的 距离 在 逐渐 缩 
短 , 在 7.2 ms 时 两 发 回转 体 发 生 追 尾 。 整 体 上 , 首 


dxiv 2 


减 。 在 入 水 撞击 的 瞬时 ,次 发 回转 体 受到 的 加 速度 
大 小 明显 小 于 首发 回转 体 ,主要 是 由 于 次 发 回转 体 
撞 水 时 ,水 层 下 方 存在 首发 回转 体 的 低压 空 泡 流 场 ， 
减弱 了 次 发 回转 体 撞 水 的 反 向 冲击 作用 。 随 着 时 间 
的 增加 ,次 发 回转 体 进入 了 首发 回转 体 的 空 泡 内 部 
运动 ,其 受到 的 阻力 急剧 减 小 ,所 以 加 速度 值 急剧 变 
小 和 速度 衰减 变 缓 。 随 后 次 发 回转 体 的 加 速度 曲线 
出 现 增 大 、 减 小 、 再 增 大 和 减 小 的 趋势 ， 这 是 由 于 随 
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Fig.11 When At=1 ms,variation of dimensionless 
penetration depth , velocity and acceleration 


of revolved bodies with time 
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图 12 At=3 ms 时 ,回转 体 无 量 纲 入 水 深度 、 
速度 、 加 速度 随时 间 的 变化 
Fig.12 When At=3 ms,variation of dimensionless 
penetration depth , velocity and acceleration 
of revolved bodies with time 
图 13 表示 的 是 入 水 时 序 为 5 ms 时 ,串联 回转 
体 无 量 纲 人 水 深度 .速度 .加 速度 随时 间 的 变化 曲 
线 。 在 1=9 ms 时 刻 , 首 发 回转 体 的 运动 速度 已 经 衰 
减 到 初速 度 的 27.6% ,此 时 次 发 回转 体 和 首发 回转 
体 还 没有 发 生 追 尾 。 在 人 水 撞击 瞬间 ,次 发 回转 体 
的 加 速度 是 小 于 首发 回转 体 的 ,与 人 水 时 序 为 3 ms 
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着 首发 回转 体 空 化 数 的 持续 增 大 及 次 发 回转 体 在 首 
发 空 泡 内 运动 的 干扰 ,引起 首发 空 泡 内 部 燕 汽 相 液 


时 的 工 况 一 样 。 在 和 人 水 阶段 , 受 首发 回转 体 运动 过 
后 留 下 的 脱落 空 泡 影响 ,次 发 回转 体 的 速度 衰减 程 


度 是 小 于 首发 回转 体 的 。 
3 一 首发 入水 


= 多 故人 水 


HID 


tms 
(a) 无 量 纲 入 水 深度 


一 一 首发 入 水 
= 次 发 入 水 


20 
0 1 2 3 4 | 6 0 8 2 
tms 
(b) 速度 
=33 
一 一 首发 入 水 
到 1 -次 发 和 水 


oa(103m。Ss) 


tms 


(c) 加 速度 
图 13 At=5 ms 时 ,回转 体 无 量 纲 入 水 深度 、 
速度 、 加 速度 随时 间 的 变化 


Fig.13 When At=5 ms,variation of dimensionless 


penetration depth , velocity and acceleration 


of revolved bodies with time 


通过 上 述 的 分 析 , 可 以 发 现 不 同 的 入 水 时 序 下 
串联 回转 体 将 出 现 不 同 特征 的 扰动 ,该 扰动 对 回转 
体 入 水 空 泡 的 发 展 和 运动 特性 产生 重要 的 影响 。 随 
着 入 水 时 序 的 增加 ,依次 出 现 两 种 运动 模式 ,分 别 为 
追尾 模式 和 弱 扰 动 模式 。 追 尾 模 式 为 次 发 回转 体 受 
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到 首发 回转 体 空 泡 流 场 的 影响 ,次 发 回转 体能 够 进 
入 到 首发 回转 体 的 空 泡 内 部 运动 ,最 终 次 发 回转 体 
将 追 上 首发 回转 体 。 弱 扰动 模式 为 首发 回转 体 的 人 
水 空 泡 演化 未 受 影响 ,而 次 发 回转 体 的 空 泡 演化 和 
运动 受到 首发 回转 体 空 泡 尾 迹 的 干扰 。 


论 


3 结 
本 研究 建立 了 模拟 串联 回转 体 和 人 水 过 程 的 三 维 
数值 模型 , 开展 了 不 同人 水 时 序 下 串联 回转 体 以 
106. 8 m/s 的 初速 度 入 水 问题 研究 ,获得 了 典型 工 况 
下 串联 回转 体 的 入水 空 泡 演化 特性 以 及 串联 回转 体 
无 量 纲 和 水 深度 .速度 .加 速度 随时 间 的 变化 。 主 要 
获得 如 下 结论 。 
志 1) 次 发 回转 体 入 水 过 程 受 到 首发 回转 体 和 水 空 
泡 滴 场 的 影响 ,在 入 水 撞击 阶段 ,次 发 回转 体 受到 的 
着 王 力 小 于 首发 回转 体 。 而 且 在 人 水 过 程 中 次 发 回 
转 体 的 速度 衰减 小 于 首发 回转 体 。 
2 ) 随 着 人 水 时 序 增 大 ,串联 回转 体 人 水 运动 模 
式 柜 次 出 现 追 尾 模 式 和 弱 扰 动 模式 。 在 追尾 模式 
中 次 发 回转 体能 够 进入 首发 回转 体形 成 的 入 水 空 
涝 内 部 运动 , 当 入 水 时 序 较 小 时 ,次 发 回转 体 不 产生 
/ 情 空 泡 ; 当 和 水 时 序 较 大 时 ,首发 次 发 回转 体 各 
自 产生 的 空 泡 会 融合 成 一 个 空 泡 。 在 弱 扰动 模式 
首发 回转 体 的 空 泡 演化 过 不 受 次 发 回转 体 入 水 
的 影响 ,而 次 发 回转 体 的 空 泡 演化 和 运动 受到 了 首 
发 问 转 体 空 泡 尾 迹 的 干扰 。 
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